Surface electronic structure and magnetic properties of doped manganites 
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The electronic structure and magnetic properties of Lai-^A^MnOa are investigated. It is assumed 
that, at the outermost layer, the environment of the Mn ions does not have cubic symmetry. The 
e g orbitals are split and the double exchange mechanism is weakened. The charge state of the Mn 
ions is modified, and the magnetic ordering of the spins tends to be antiferromagnetic. The surface 
magnetization and the dependence of the transport properties through the resulting surface barrier 
on applied magnetic field and temperature is analyzed. 
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Doped manganites show many unusual features, the 
most striking being the colossal magnetoresistance in 
the fully ferromagnetic phase [0J|]. Extensive research 
shows that the transport properties, and the magne- 
toresistance in particular, are significantly modified at 
artificially created barriers p|-|TcH or in ceramic materi- 
als P~q1 . The magnetoresistance, as function of tem- 
perature, drops more rapidly than in the bulk. It has 
larger values at low fields, and persists at large fields, 
unlike in the bulk case. The relevance of the interfaces in 
perovskite manganites can also be inferred by comparing 
with transport in related materials which exhibit colossal 
magnetoresistance ]16[ |. These properties have been as- 
cribed to changes in the interface structure, although the 
origin of these modifications, and the resulting structure 
are not known. 

In the present paper, we analyze the simplest, and 
most common, modification with respect to the bulk that 
a surface may show: the loss of cubic symmetry around 
the Mn ions. It is well known that Lai-^A^MnOa shows 
a transition from a tetragonal (or orthorhombic) struc- 
ture to a cubic one as the value of the doping x is in- 
creased. The systems with the highest Curie tempera- 
ture have x ~ | and are in the cubic phase. This implies 
that the two e g orbitals of the Mn ion are degenerate and 
contribute to the conduction band. In this situation, the 
double exchange mechanism is enhanced. 

The cubic symmetry is lost at the surface. If the last 
layer is oxygen deficient, the resulting splitting between 
the e g orbitals can be larger than typical Jahn- Teller 
splittings in Lai-^A^MnOs with small values of x. When 
one of the e g orbitals moves away from the Fermi level, 
the double exchange mechanism is weakened, and direct 
antiferromagnetic couplings between the core S = | spins 
can prevail. Moreover, the reduction in electronic kinetic 
energy can also lead to charge transfer between the sur- 
face layers and the bulk, contributing to the formation of 
a surface dipole. All these effects can be modified by sur- 
face spin waves, which, in turn, depend on temperature 
and external magnetic fields. A significant dependence 
of the metal-insulator transition temperature as function 
of oxygen pressure in thin films is found in [fl7| . 

In order to investigate these features, we start from 



a tight binding Hamiltonian, using the two e g orbitals, 
d x 2_ y 2 and d 3z 2_ r 2, which we designate x and z respec- 
tively. Hopping between them takes place through vir- 
tual jumps to the intermediate O ions. Fixing the ori- 
entations of the e g orbitals to the frame of reference of 
the lattice, we obtain for the z— direction, t zz = t and 
t X z = txx = 0. For the directions in the x — y plane 
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and t xx — |t, where the two signs 
correspond to the x and y directions, and t is the 



hopping generated from the (dpa) matrix 
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element between a d orbital in a given Mn ion, and ap or- 
bital in a neighboring O ion JlJ,|l^l . At each site there is 
also a spin, from the three ti g orbitals, which we treat as 
classical, and parametrize in terms of the Euler angles 
and (f>. We assume that the Hund's coupling between the 
e g electrons and this spin is much larger than other scales 
in the model. Then, the hopping elements depend on the 
orientation of the core spins, and the actual hopping is 
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where the value of the hopping has been estimated t ~ 
Q.l-0.3eV. 

We study a cubic lattice with periodic boundary condi- 
tions in the x and y directions, and open boundary con- 
ditions along the z direction, where two (001) surfaces 
terminate the lattice. We use slabs with a thickness of 
20 atoms. We have checked that, for this size, the bulk 
properties are recovered at the center of the slab. The 
environment of the outermost Mn ions is deficient in oxy- 
gen. The oxygen octahedra which surround the Mn ions 
are incomplete. The absence of the negatively charged 
2 ~ ions leads to a downward shift of the e g levels with 
respect to the values in the bulk. This shift is larger 
for the d^ z 2_ r 2 orbital, which points towards the surface. 
The d x 2_ y 2 orbital is more localized around the ion, and 
is less sensitive to the change in the environment. In or- 
der to keep the number of free parameters in the model to 
a minimum, we leave the d x 2_ y 2 level unchanged with re- 
spect to the bulk, while the d% z 2_ r 2 is shifted downwards 
by an amount A. The value of A should be compa- 
rable, or larger, than the observed Jahn- Teller splitting 
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in LaMn03, which, in turn, is larger than the e„ band- 
width. Reasonable values of A are ~ 0.5 — 1.5eV |2J]. The 
e g levels at all other layers remain unchanged, except for 
electrostatic effects (see below). 

This shift of the surface d 3z 2_ r 2 orbitals leads to charge 
transfer between the bulk and the surface. We treat the 
induced electrostatic effects within the Hartree approx- 
imation. This gives rise to an additional shift in the 
electronic levels which is determined by solving selfcon- 
sistcntly the Schrodinger and the Poisson equations. This 
electrostatic shift is equal for the two e g orbitals at each 
layer. Screening from other levels is described in terms 
of a dielectric constant of value e = 5eo- Selfconsistency 
is imposed on the e g levels in all layers in the system. 
We do not consider the possibility of orbital ordering in- 
duced by a Hubbard U between electrons at the two e g 
orbitals We analyze underdoped materials, x < 0.5, 
where no unusual magnetic textures are expected p2| . 
The Hartree approximation acts to suppress charge fluc- 
tuations, although it does not split the e g bands in the 
way a Hubbard term does. This effect should be less im- 
portant at the surface, due to the crystal field splitting 
A introduced before. 

Typical results, for A = lOt, hole concentration of 
x=0.3 and a ferromagnetic configuration of all spins are 
shown in fig. jj] . Charge neutrality corresponds to a total 
occupancy of the e g orbitals of 0.7. 



ues of A/t do not alter these results. For intermediate 
values of A/t, the surface electrostatic barrier depends 
significantly on the magnetic configuration. The surface 
d x 2_ y 2 orbitals are practically empty, and the d 3z 2_ r 2 is 
occuppied for A > 3t. 

The suppression of the double exchange ferromagnetic 
coupling at the surface leads to an enhancement of the 
superexchange interaction among the Mn ions. To in- 
vestigate further this effect, we study magnetic configu- 
rations where the core spins of the surface Mn ions are 
allowed to rotate, as shown in fig. ||. The canting an- 
gle 9 is used as a variational parameter. 9 = implies 
a perfect ferromagnetic order at the surface. 9 = ir/2 
gives rise to an antiferromagnetic alignment of the sur- 
face spins, at right angles with the bulk magnetization. 
To take into account the superexchange interaction, we 
introduce an antiferromagnetic coupling between the Mn 
core spins. Substracting a trivial constant, the energy 
per surface ion, due to this coupling, is: 



E AF = J ± cos(0) + 2 Jy cos(26») (2) 
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FIG. 1. Charge in the e g orbitals as a function of the dis- 
tance from the surface for A = lOt. On the surface, ul, = 1, 
the x orbitals are almost empty, due to the shift A. 

The large splitting between the e g orbitals at the sur- 
face leads to a significant charge transfer to the outermost 
layer. The charge distribution reaches the bulk values at 
the third layer, in agreement with the expected short 
screening length of the metal. The surface d^ z 2_ r 2 or- 
bitals are almost full, while the d x 2_ y 2 are empty, so that 
the charge state of the ion is Mn 3+ . Hence, for the value 
of A used in fig. Q, the double exchange mechanism is 
almost completely suppressed at the surface. Higher val- 



where 9 is the angle shown in fig. &, J± is the antifer- 
romagnetic coupling between a Mn spin at the surface 
and one in the next layer, and Jy is the coupling between 
spins at the surface layer. From symmetry considera- 
tions, J|| = 4J ± Q. 
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FIG. 2. Magnetic structure considered in the text. Only 
the spins at the surface layer are rotated, in the way shown. 



The total energy of the system is the sum of the kinetic, 
Hartree and magnetic energies. By minimizing the total 
energy as a function of 9 we obtain the canting angle as 
a function of the antiferromagnetic coupling J± . 

Fig. H shows the calculated phase diagram as func- 
tion of the splitting at the surface between the e g levels, 
A, and the direct antiferromagnetic coupling between the 
core spins J± , for two different values of the hole concen- 
tration. 
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FIG. 3. Magnetic phase diagram of the surface of a doped 
manganite. The lines separate the fully ferromagnetic and 
the canted regions. The dots indicate that the order is 90% 
antiferromagnetic, namely, 6 — 81°. Full dots are results for 
hole concentration of x = 0.3. Open dots are for x = 0.2. 
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FIG. 4. Magnetization of the last layer as function of tem- 
perature for a cluster with 20 layers and open boundary con- 
ditions. For comparison, the results for a cluster of the same 
size and periodic boundary conditions (no surfaces) are also 
shown. 



Realistic values of Jj_ are of the order of 0.02i j|24f] 
and from fig. 0, we see that, for this value, the value of 
9 is close to that of perfect surface antiferromagnetism, 
7r/2. Our approach tends to underestimate this tendency 
towards antiferromagnetism, as we do not allow to relax 
the spins in the layers deeper into the surface. 



The changes in the magnetic surface structure also lead 
to modifications in the spin stiffness at the surface, which 
is weaker than in the ferromagnetic bulk. It is straight- 
forward to show that weaker couplings lead to the for- 
mation of bands of surface magnons pq ]. We find such 
bands both in the canted and in the ferromagnetic phases 
shown in fig. ||. The bandwidth is narrower than the 
bulk magnon band. Thermal excitation of these modes 
leads to a decay of the surface magnetization as func- 
tion of temperature which is faster than in the bulk. In 
order to estimate this effects, we have calculated, using 
Monte Carlo techniques [pGfl , the surface magnetization 
of a cluster of classical spins. In this model, the bulk dou- 
ble exchange mechanism is described by an effective fer- 
romagnetic Heisenberg coupling J, and the surface spins 
interact with an antiferromagnetic coupling J\ . The out- 
ermost spins interact with the spins in the second lay- 
ers via a ferromagnetic coupling J 2- From the results 
reported previously, we estimate that J\ is very small, 
r-j — J] 100 and J2 takes a value between J/2 and J/6. In 
fig. jj] we plot the surface magnetization as a function 
of the temperature for different values of J2 . The results 
obtained do not depend strongly on the value of J2 and 
they are in reasonable agreement with the experimental 
results reported in M. 



The combination of energy shifts and changes in the 
magnetic couplings reduces the transmission of the sur- 
face, even in the absence of other, extrinsic, barriers. In 
order to estimate this effect, we have calculated the re- 
sistance between two perfect double exchange ferromag- 
netic metals separated by a layer with the amount of 
charge and the magnetic structure obtained in the pre- 
vious calculations. Since only the d 3z 2_ r 2 orbitals con- 
tribute to the transport through the interface, we have 
simplified the model, and only a single orbital per site 
is kept. The calculations were done using the method 
described in [^,^|. Two effects contribute to the resis- 
tivity, the shift in energy of the interface orbital and the 
difference in spin orientation between the atoms at the 
bulk and at the interface. The conductance is sharply 
reduced when the surface layer is antiferromagnetic, as 
the double exchange mechanism is suppressed. We find 
that for hole concentrations of x=0.3 the presence of the 
interface increases the resistance of the system in a factor 
bigger than 10. 

A magnetic field will reduce the antiferromagnetism 
at the surface, and, at high fields, the surface spins are 
aligned parallel to the bulk. We estimate this effect by 
adding a magnetic field to the model and finding the 
magnetic structure which minimizes the energy. The re- 
sulting conductance, for A = 3t and J± = 0.025i, is 
plotted in fig. ||. A very high field ~ 60T is required to 
saturate the magnetoresistance. Note that the low field 
( j IT ) magnetoresistance is probably due to the align- 
ment of the polarization of the bulk electrodes, which can 
be understood within conventional models 1 30 3^]. Our 
results for the high field dependence are consistent with 
the available experimental data data |^10[. 
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FIG. 5. Inverse conductance, normalized to the conduc- 
tance of a perfect ferromagnetic system (8 = and A = 0), 
of a magnetic layer, as function of applied field. The arrow 
indicates the saturation limit at high fields. The size of the 
system is 20 x 20 x 20. The parameters used are described in 
the text. 

The magnetic excitations of the surface layer also mod- 
ify the magnetoresistance at finite temperatures |32] Q . 
Spin flip scattering due to thermally excited magnons 
leads to a suppression of the magnetoresistance at tem- 
peratures below the bulk Curie temperature. The depen- 
dence of the magnetoresistance on temperature should 
be similar to that of the surface magnetization, shown in 

In conclusion, we have shown that the lack of cubic 
symmetry at surfaces, combined with the double ex- 
change mechanism, leads to significant changes in the 
magnetic and transport properties of doped manganites. 
Charge is transferred from the bulk to the surface layers, 
and an antiferromagnetic ordering of the surface spins is 
favored. This structure is consistent with the observed 
reduction of the magnetoresistance in junctions and ce- 
ramic materials at temperatures below the bulk Curie 
temperature, with the high field dependence of the mag- 
netoresistance and with the experimental results on sur- 
face magnetization. 
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